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Abstract The Bel- Robinson tensor is analyzed as a linear map on the space 
of the traceless symmetric tensors. This study leads to an algebraic classifi- 
cation that refines the usual Petrov-Bel classification of the Weyl tensor. The 
5_i ' new classes correspond to degenerate type I space-times which have already 

OJQ. been introduced in literature from another point of view. The Petrov-Bel 

types and the additional ones are intrinsically characterized in terms of the 
sole Bel-Robinson tensor, and an algorithm is proposed that enables the dif- 
ferent classes to be distinguished. Results are presented that solve the prob- 
lem of obtaining the Weyl tensor from the Bel-Robinson tensor in regular 
OO . cases. 



Keywords Bel-Robinson tensor • Gravitational superenergy • Petrov-Bel 
classification 



1 Introduction 



With the aim of defining intrinsic states of gravitational radiation Bel Q] [2] 
[3] introduced a rank 4 tensor which plays an analogous role for gravitational 
field to that played by the energy tensor for electromagnetism. This super- 
energy Bel tensor is quadratic in the Riemann tensor and, in the vacuum 
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case, is divergence-free. In this case the Bel tensor coincides with the super- 
energy Bel-Robinson tensor, built with the same expression by replacing the 
Riemann tensor with the Weyl tensor. 

Both, the Bel and the Bel-Robinson tensors have interesting properties 
that, relative to every observer, allow the definition of a non negative super- 
energy density and a super-energy Poynting vector. In the last decade the 
interest in the super-energy tensors has been on the increase, and a lot of 
works are devoted to analyzing their properties and to studying their gener- 
alization to any dimension and to any physical field [I] (see references therein 
for an exhaustive bibliography on this subject). In a recent work [5], where 
the dynamical laws of super-energy are accurately analyzed, up-to-date ref- 
erences on the Bel and Bel-Robinson tensors can be found. 

The algebraic properties of the Bel-Robinson tensor (BR tensor) and its 
close relationship with the principal null directions of the Weyl tensor were 
studied by Debever early on [6] [7], and a spinorial approach can be found 
in [8|. But the intrinsic algebraic characterization of a BR tensor was not 
obtained until the recent work by Bergqvist and Lankien [9] . They give the 
conditions on the BR tensor playing a similar role to that played by the 
algebraic Rainich [Hj] conditions for the electromagnetic field. 

Nevertheless, some gaps remain in the algebraic comprehension of the BR 
tensor. The analogy with the electromagnetic field can help us to understand 
them. Let us consider an energy tensor T em satisfying the algebraic Rainich 
conditions: 

1 



T em=X 2 9, X= 2^J^ T em- (1) 

For a T em satisfying fl} wc know (see for example [H]): 

(i) Its algebraic classification: when x / 0, the Segre type of T em is 
[(11)(11)] and it is associated with a non null electromagnetic field F; when 
X = 0, the Segre type of T em is [(31)] and it is associated with a null electro- 
magnetic field F. 

(ii) Its canonical expression in terms of its intrinsic elements, namely, 
invariant subspaces and scalars: for a non null field, T takes the expression 
T em = — xH, where 77 is a 2+2 structure tensor that determines the two 
principal planes and the two null principal directions of the electromagnetic 
field; for a null field, T takes the expression T em = I Cg) I, where I is the 
fundamental vector of the electromagnetic field. 

(iii) It determines the electromagnetic field F up to a duality rotation, 
and the explicit expression of F in terms of T em is also known |llj . 

As commented above, Bergqvist and Lankien |9J have found the necessary 
and sufficient conditions for a rank 4 tensor T to be the BR tensor associated 
with a Weyl-like tensor W . These conditions play for the BR tensor the same 
role as the Rainich conditions (p} play for the electromagnetic field. On the 
other hand, the three points satated above for the electromagnetic field could 
be similarly established for the BR tensor: 

(I) Algebraic classification of a BR tensor T: the Petrov-Bel classification 
can be obtained by studying the Weyl tensor as a linear map [T2], [3]. The 
symmetries of the BR tensor T allow us to consider and analyze it as a linear 
map on the nine dimensional space of the traceless symmetric tensors. What 



algebraic classification follows on from this study? What relationship exists 
between these BR classes and the Petrov-Bel types of the Weyl tensor? 

(II) Canonical form of the BR tensor T in terms of its invariant spaces 
and scalars: for every algebraic type of the BR tensor, the eigenvalues and 
eigenvectors should be analyzed, as well as the canonical expression of T in 
terms of them. What relationship exists between the spaces and scalar the 
invariants of both, the BR tensor and the Weyl tensor? 

(III) Expression of the Weyl tensor in terms of the BR tensor: it is known 
that the BR tensor T determines the Weyl tensor W up to a duality rotation. 
But the explicit expression of W in terms of T has not been established. 

The goal of this work is to tackle the algebraic problems of the BR tensor 
stated in the three points above. Here we solve the first one and give prelim- 
inary results on the points (II) and (III) which will be fully solved in another 
work in progress [13] . 

Some basic properties of the BR tensor T considered as a linear map 
are presented in section [2] of this work. In section [3] we define the main BR 
invariant scalars and study their relationship with the scalar invariants of 
the Weyl tensor. This result is a contribution to point (II) above. We also 
study the powers of T and obtain a family of identities on the BR tensor that 
generalize an already known equality. Finally, we determine the characteristic 
polynomial of T . 

In section 0] we show that the main BR scalar invariants allow different 
classes of Weyl tensors to be distinguished. On one hand, they label the three 
groups of Petrov-Bel types that can be discriminated by the Weyl eigenval- 
ues, namely, types III, N and O with a triple Weyl eigenvalue, types II and 
D with a double eigenvalue and type I with three different eigenvalues. On 
the other hand, the main BR scalars also distinguish between five classes of 
algebraically general Weyl tensors. These classes refine the Petrov-Bel clas- 
sification and were introduced by Mcintosh and Arianrhod [14] . They corre- 
spond to particular configurations of the four null principal directions: either 
they span a 3-plane or they define a frame with permutability properties. A 
detailed analysis of these 'degenerate' type I classes can be found in [15] , 

In section [5] we consider the classification of the BR tensor taking into 
account its eigenvalue multiplicity. We show that seven classes appear from 
this point of view: one with 9 different eigenvalues, two with 6 different 
eigenvalues, three with 3 different eigenvalues and one with a eigenvalue of 
multiplicity 9. We characterize every class in terms of the main BR invariant 
scalars and show that they precisely correspond to the classes of the Weyl 
tensor considered in the previous section. An arrow diagram that helps us to 
visualize the different classes and their degenerations is offered at the end of 
the section. 

The study presented in sections Q] and \5\ shows that the main BR scalars 
do not distinguish between Petrov-Bel types III, N and O, and between types 
II and D. Nevertheless we prove in section [7] that all these Petrov-Bel types 
match with different classes of the BR tensor and that they can be discrimi- 
nated by using their minimal polynomial. At the end of this last section we 
present an algorithm that summarizes the classification of the BR tensor and 
the intrinsic characterization of every class. 



In studying the minimal polynomial in last section we use a expression 
that relates the n-powers of the BR and Weyl tensors. This expression is 
proved in section [6] and it is valid except for some degenerate cases. For the 
case n = 3 this expression allows us to obtain, in the same section, the Weyl 
tensor in terms of the BR tensor when none of the main BR scalar invariants 
vanish. This result partially solves point (III) above. Elsewhere [13] we obtain 
the Weyl tensor in terms of the BR tensor in the cases that are not considered 
here. 

With the intention of making the mathematical expressions simpler we 
use, when possible, a global notation that is accurately explained in four 
appendixes at the end of this work. This also allow us to simplify the calcu- 
lations by using some properties that are also stated in these appendixes. 

It is worth outlining that, in spite of the fact that the BR tensor loses 
part of the information of the Weyl tensor (a duality rotation), its study as 
a linear map leads to ten different classes, a richer classification than the six 
Petrov-Bel types of the Weyl tensor as a linear map. On the other hand, note 
that from our study follows a characterization of every Petrov-Bel type in 
terms of the BR tensor, a result that only was known for types N and O [16] 



2 The Bel-Robinson tensor: towards its algebraic classification 

We shall note g the space-time metric with signature { — ,+,+,+} and we 
shall write W to indicate the Weyl tensor. The conventions of signs are those 
of the book by Stephani et al. [18] . The Bel- Robinson tensor (BR tensor) T 
is given in terms of the Weyl tensor as: 

T a ^„ = \ (W a p /W„ P ua + *W a V * WW) » (2) 

where * denotes the Hodge dual operator (see Appendix 1519). Note that 
expression (|2|) coincides with that given in [8] [17] and it differs in a factor 
from the original one by Bel [3]. 



2.1 The Weyl tensor as an endomorphism 

A self-dual 2-form is a complex 2-form T such that *J- = \T . We can 
associate biunivocally with every real 2-form F the self-dual 2-form T = 
A={F — i * F). In short, here we refer to a self-dual 2-form as a SD bivector. 
The endowed metric on the 3-dimensional complex space of the SD bivectors 
is Q = h{G — i rf), r\ being the metric volume element of the space-time, and 
G being the metric on the space of 2-forms, G = \g Kg. Here A denotes the 
double-forms exterior product (see Appendix |A"15). 

The algebraic classification of the Weyl tensor W can be obtained [12], 
[3] by studying the traceless linear map defined by the self-dual (SD) Weyl 



tensor VV = \ (W — i * W) on the SD bivectors space. This SD-endomorphism 
(see notation in Appendix [B]) has associated the complex scalar invariants: 

a = TrW 2 , 6 = TrW 3 . (3) 

In terms of them, the characteristic equation reads 

1 1, 
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6 = 0. (4) 



Then, the Petrov-Bel classification follows taking into account both the eigen- 
value multiplicity and the degree of the minimal polynomial. The algebraically 
regular case (type I) occurs when the characteristic equation (gj) admits three 
different roots. If there is a double root p = — - and a simple one — 2p, the 
minimal polynomial distinguishes between types D and II. Finally, if all the 
roots are equal, and so zero, the Weyl tensor is type O, N or III, depending 
on the minimal polynomial. 



2.2 The Bel-Robinson tensor as an endomorphism 

The expression ([2]) of the BR tensor may be written in terms of the SD Weyl 
tensor as: 

T a ^ v = W a p a W lipV( r , (5) 

where ~ stands for complex conjugate. Taking into account the o-product 
defined in Appendix [D] we have that (O becomes 

T = WoW. (6) 

It is easy to show that, for two arbitrary SD-endomorphisms X and y, the 
tensor Q = X o y satisfies 

QafifSu = Qfj.afli' = Qflvotfj, j Q a f3u = " ■ {') 

As a consequence of these properties, if S is a trace-less symmetric tensor 
(TLS tensor), then Q(S) a p = Q a ^ v S^ is trace-less and symmetric. Thus, 
we can see Q as a (complex) symmetric linear map on the space of the TLS 
tensors, that is to say, Q is a TLS- endomorphism (see notation in Appendix 

EJ). 

Besides, the BR tensor as given by (JSJ) has two important additional 
properties. Firstly, it is real and, secondly, as a consequence of the traceless 
property of the SD Weyl tensor, it is a fully symmetric traceless tensor (see 
Appendix [DJ 7). Then we have: 

Lemma 1 The Bel-Robinson tensor T is a real traceless TLS- endomorphism, 
that is, a traceless symmetric endomorphism on the 9-dimensional space of 
the traceless symmetric tensors. 



As explained in subsection above, the Petrov-Bel classification follows on 
by studying the Weyl tensor as a SD-endomorphism. Similarly, we want to 
obtain the algebraic classification of the BR tensor by analyzing it as a TLS- 
endomorphism. We start our study by writing the BR tensor in terms of 
another TLS-endomorphism. If we define 

n = w°g, (8) 

and we take into account ^ and that T is real, we obtain: 

T= f2»J] = ~n»n, (9) 

where • denotes the composition of TLS-endomorphisms (see Appendix [Cj . 
The tensor fl is not but the restriction of the Weyl tensor acting (with the 
indices 2 and 4) on the TLS tensor space. Indeed, for a trace-less symmetric 
tensor S, we have 

Q^P^ S^ v = w a / V S"". (10) 

On the other hand, it is easy to show that the tensorial components of fl are 

^l a f3^ = -W a (p,iy)(3 7 (H) 

where ( ) denotes symmetrization. 

The tensor fl inherits some properties of the SD Weyl tensor. Evidently, 
it is a TSL-endomorphism (it satisfies (|Tj) ) and, moreover, it is traceless and 
satisfies a SD-like property: 

^ Q /iQiy = 0' ^a-tpsG fiu = 0. (12) 

Conversely, if fl is a (complex) TLS-endomorphism that satisfies (fl"2|) we 
can reverse the expression (|5J| and obtain the traceless SD-endomorphism 
(see property of Appendix |Pl5): 

2 4 

Wapftv = - fl a ~tI38 Q 1 fiv = o^a[H/3 ' (^) 

where [ ] denotes antisymmctrization. Thus we have the following. 

Proposition 1 There is a one-to-one correspondence between traceless SD- 
endomorphisms W and the (complex) TLS-endomorphism fl which satisfy 
(fi~2|) . This bijection is given by ([8]) and (fl~3|) . 



2.3 Can the Weyl tensor be obtained from the Bel-Robinson tensor? 

The BR tensor can be obtained from the Weyl tensor by means of the 
quadratic expression ([2|), but it is known that this expression is invariant 
under duality rotation of the Weyl tensor. For the SD Weyl tensor a dual- 
ity rotation takes the form e W. Then, the above quoted invariance follows 
trivially from the expression ([6]) because, 



WoW= [e ie W}o[e- ie W]. 



Thus, the accurate question we should pose is: can the Weyl tensor be 
determined, up to a duality rotation, from the BR tensor? or, to be more 
precise, is there an explicit algorithm to obtain it? These questions will be 
analyzed and solved elsewhere 13J for an arbitrary BR tensor. In section[S]we 
study this question for regular cases. In this study the TLS endomorphism 
fl plays an important role. 



3 Powers, scalar invariants and characteristic polynomial of the 
Bel-Robinson tensor 

In order to tackle the algebraic classification of the BR tensor T as a TLS- 
endomorphism we need to know its characteristic polynomial and its scalar 
invariants. In this section we obtain them from the powers of T. 

The powers of the BR tensor can be obtained as T n = T n ~ 1 • T. From 
© and as T is real, the fl defined in ([5]) satisfies Q» i? = flu f2. Then, from 
(J9j) we obtain 

Lemma 2 For n ^ 1 it holds T n = Q n • il , where [2 is the endomorphism 
defined as ((SJ) from the Weyl tensor W. 

On the hand we can obtain the TLS-endomorphism Q n = Q • i?™ -1 in 
terms of the SD-endomorphism W™ = W oYV n ~ l . Indeed, from the definition 
© and putting, in expression j38|) of Appendix [Dj8, V = W, X = W"" 1 
and Z = y = Q , we obtain by induction: 

Lemma 3 The TLS-endomorphism fl given by ((5J satisfies 

n n = w n o~g, n >i. (14) 

The expression for the powers of the BR tensor and the traces of these 
powers in terms of the SD Weyl tensor can also be obtained by induction by 
taking, in expression ([!§ of Appendix |Ej]8, V = Z = W, X = y = W n_1 , 
and taking into account the property of Appendix [Dl2. Thus, we have: 

Proposition 2 The powers of the BR tensor and their traces can be com- 
puted in terms of the SD Weyl tensor as 

r = W"oW", trT" = |trW™| 2 . (15) 

From this proposition, and by applying the property of Appendix [Dj2, we 
obtain the following result easily: 

Corollary 1 The powers T" of the BR tensor satisfy: 

(T n ) X a , P = \TrT n g ap . (16) 

For n = 2 this expression states that T a x^Tp Xfn " = jT pX ^T pXf11 - 'g a0 , an 
identity already known [6] [8]. Thus (| 16[) generalizes, for an arbitrary n, this 
identity. 

As a consequence of proposition [2j the traces of the powers of T are non 
negative and we can compute them in terms of the traces of the SD Weyl 



tensor. But we also know that W satisfies the characteristic equation ((4J, 
that is: 

W 3 = ^aW+^bG, (17) 

where a and b are the main Weyl scalar invariants given in ([3]). 

From (fT"7)) . the traces of the powers of W™ for n > 3 can be computed in 
terms of a and b. After that, they can be used to compute the traces of T n 
by applying proposition O More precisely, if we denote 

(18) 



a = -\a\, 


/3 = 


i> 


M = 


38. (■*■ + **), 


tave that the traces of T n 


are j 


jiven by 




TrT = 
TrT 3 = 9/3 2 
TrT 5 = 25a 2 /3 2 
TrT 7 = 49a 4 /3 2 
TrT 9 = 9/3 2 [/3 4 + 9a 


6 + 3/i] . 


TrT 2 = 
TrT 4 = 
TrT 6 = 
TrT 8 = 


:4q 2 

:4q 4 

: 4q 6 + 9/3 4 + 6^ 

= 4q 8 + 64a 2 /3 4 + 12a 2 /* 



(19) 



Note that the Weyl tensor defines four real scalar invariants that can be 
grouped in the complex ones a, b defined in © from the SD Weyl tensor. 
We have remarked in subsection 12.31 that a duality rotation is lost when 
constructing the BR tensor. So just three scalars survive in T: we can see in 
fl!9p that two of them, the modulus of a and b determine the traces of T 2 
and T 3 . The third one, a?b 2 +a 3 6 2 , does not appear until TrT 6 is computed. 
This fact and expressions (fT9|) justify defining the following main BR scalar 
invariants: 



a 



IVTiT 2 , [3=IVTyT\ u= iftrT 6 -4a 6 -9/3 4 V (20) 

2 3 6 V / 



Then, the traces of the powers of the BR tensor depend on the BR main 
scalar invariants as ([19]) . Moreover, these expressions can be used to obtain 
the characteristic equation of the BR tensor and we achieve the following. 

Proposition 3 The BR tensor satisfies the characteristic equation 

x 9 - 2a 2 x 7 - 3fJ 2 x 6 + a 4 x 5 + a 2 (3 2 x i + (3/5 4 - fi)x 3 + a 2 (3 i x - (3 6 = (21) 

where a, (3 and \x are the main BR scalar invariants (|20|) . These invariants 
determine the traces Tr T n as (TH 



4 Labeling some Weyl types with the main BR scalar invariants 

In the next section we study the algebraic types of the BR tensor as TLS- 
endomorphism. In order to understand the correspondence between these 
new BR types and the already known Weyl types, we show in this section 
that some of these Weyl classes can be characterized in terms of the main 
BR scalar invariants. 



When analyzing the Weyl tensor as a SD-endomorphism, the main Weyl 
invariant scalars distinguish the multiplicity of the Weyl eigenvalues: the 
algebraically general type I has 3 simple eigenvalues, types II and D have 
a simple and a double one, and types III, N and O have a triple vanishing 
eigenvalue. Then, from the characteristic equation (j4|) we obtain [18] |19j : 

Lemma 4 Let a — TrW 2 , b = TrW 3 be the main Weyl scalar invariants. 
Then, the Weyl tensor is 

i. Petrov-Bel type I iff 6b 2 — a 3 ^= 0. 

ii. Petrov-Bel types D or II iff 6b 2 = a 3 ^ 0. 

Hi. Petrov-Bel types 0, N or III iff a = = b. 

Note that in order to distinguish between types D and II as well as between 
types N or III, the minimal polynomial is necessary [19] . But, at the moment, 
we are only interested in scalar conditions. 

On the other hand, some classes of algebraically general space-times have 
been considered in literature. A refinement of the Petrov-Bel classification is 
based on the Weyl scalar invariant [H] : 



M= ¥ -6. (22) 

The invariant M is not defined if b = as occurs in types III and N and 
in algebraically general space-times with a vanishing eigenvalue. But we can 
extend its validity to these cases if we put M = oo whenever b = ^ a, and 
M = 0iia = = b. 

This way, it is known [14] [20] that the cases where M is real positive or 
infinity, which are called IM + or IM°° respectively, correspond to the case of 
the Weyl eigenvalues having a real ratio, ^el (the case M — oo means 
that the Weyl tensor has a vanishing eigenvalue) . An equivalent condition in 
terms of Debever null directions has also been obtained [Tl] [2D]: M is real 
positive or infinity iff the four Debever null directions span a 3-planc. 

On the other hand, the case where M is real negative IM~ has also been 
considered. It corresponds to the property of two of the SD eigenvalues having 
the same modulus (\pi\ = \pj\ for some i j£ j). Moreover, the case M = — 6 
(or a = 0) corresponds to all the eigenvalues having the same modulus. These 
conditions have also been interpreted in terms of permutability properties of 
the null Debever directions [T5] , 

A detailed analysis of these 'degenerate' algebraically general classes can 
be found in [15] . where they are also interpreted as the space-times where 
the electric and magnetic parts of the Weyl tensor are aligned for a (non 
necessary time-like) direction. These classes also contain the purely electric 
and purely magnetic space-times which have been accurately studied in [21] . 

The scalar M is homogeneous in the Weyl tensor and thus invariant by 
duality rotation. Consequently, it should depend on the main BR scalars. 
Indeed, if we take into account the relations ([18]) between the main Weyl 
scalars a, b and the main BR scalars a, f3 and fi, we obtain: 

81/? 8 (im[Af]) 2 = 16(4a 6 /3 4 -/^ 2 ), 9 /3 4 Re[M] =4^-54/3 4 . 
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The first expression shows that M is real if, and only if, p 2 — 4a 6 /3 4 = 0. 
After that, if M is real, the second expression gives its sign in terms of the 
main BR scalars. More precisely, we have: 

Lemma 5 The invariant scalar M is real, if and only if /J? — ^a & (3 A — 0. If 
this condition holds, the Weyl tensor is: 

i. Type IM + (M real positive) iff 2p - 3 3 /3 4 > 0. 

ii. Type IM°° (M = oo) iff [3 = ^ a. 

Hi. Type IM~ (M real negative) iff 2p — 3 3 /? 4 < 

iv. Type /Ml- 6 ] (M = -6) iff a = ^ (3. 

v. Algebraically special (M = 0) iff 2p — 3 3 /3 4 = 

From now on we denote I r the class of non 'degenerate' type I Weyl 
tensors, that is to say, those with non real invariant M. Then, from lemmas 
2] and [5] and relations (JT5J) we obtain the following characterization of some 
classes of the Weyl tensor in terms of the main BR scalars: 

Proposition 4 If a, (3 and p are the main BR invariant scalars, then the 
Weyl tensor is: 

i. Type 0, N or III iff a = j3 = 0. 

ii. Type D or II iff p 2 - 4a 6 /? 4 = 0, 2 2 a 3 = 3 3 /? 2 ^ 0. 

in. Type IM+ iff y? - 4a 6 /3 4 = 0, 2fj, - 3 3 /? 4 > 0. 

iv. Type IM°° iff a^0 = /3. 

v. Type IM~ iff (J 2 - 4:a 6 f3 4 = 0, 2/i - 3 3 /3 4 < 0. 

vi. Type JMl-6] iff (3^Q = a . 

vn. Type I r iff fi 2 - Aa 6 (3 A ^ 0. 



5 Classifying the Bel-Robinson tensor 

In this section we study the BR tensor as a TLS-endomorphism taking into 
account the eigenvalue multiplicity. 

As a consequence of ©, fi and W satisfy the same minimal polynomial, 
that is, both have the same eigenvalues pi, i = 1,2,3. On the other hand, 
from expression ([9]), Q and fl are symmetric endomorphisms that commute 
and, consequently, they have the same eigenvectors. These eigenvectors and 
their associated eigenvalues should be either real or pairs of complex conju- 
gates and then T will have the same eigenvectors with associated eigenvalues 
pipj. Thus, we obtain the BR eigenvalues without solving the characteristic 
equation (|2"Tj) . 



Proposition 5 The BR tensor has 3 real eigenvalues U, and 3 pairs of com- 
plex conjugate eigenvalues t%, fi. They depend on the Weyl eigenvalues pi as 

U = \pi\ ; n=pjpk, (i,j,k) a pair permutation of (1,2,3) . (23) 

The BR tensor has three main scalar invariants. In accordance with this 
fact, the BR complex eigenvalues should depend on the three real ones. In- 
deed, a straightforward calculation shows 



11 



Lemma 6 In terms of the real BR eigenvalues t% the complex ones Tk take 
the expression: 

Tk=Pk + iq, 2pk=tk-U-tj , (i,j,k^), q 2 = P1P2+P2P3+P3P1 ■ (24) 

Generically, the BR eigenvalues are nine different scalars. Now we analyze 
the admitted degenerations. If U = tj 7^ tk, then (f24f implies pi — pj, that is, 
Tj = Tj and fj = fj and, consequently, T has 3 double eigenvalues. Moreover, 
expression (|23j) implies that the Weyl tensor has two eigenvalues with the 
same modulus, \pi\ — \p,j\, a case corresponding to a Weyl tensor of type 
IM~ . If we avoid a higher degeneration by considering tk 7^ U = tj , then 
M ^ -6. 

If one of the complex eigenvalues Tj is real, we have Tj = fj. Then (|24|) im- 
plies 5 = and consequently, tj, = fj, for every k — 1,2,3. Thus T has, again, 
3 double eigenvalues. Now, from (123| . the ratio of the Weyl eigenvalues — is 
real for every pair (i,j) because pipj is real. If we avoid higher degenerations 
we must remove the case of a vanishing Weyl eigenvalue (M 7^ 00) and then 
the Weyl tensor is type IM + . 

Thus, the case of nine different eigenvalues corresponds to a regular type 
I r Weyl tensor. From all these considerations and taking into account propo- 
sition |4] we have the followings. 

Proposition 6 A BR tensor T has nine different eigenvalues (3 real and 
3 pairs of complex conjugate) iff p? — 4a e p 4 7^ 0, that is to say, the Weyl 
tensor is type I r . 

Proposition 7 If p 2 — 4a 6 /3 4 = 0, the BR tensor T has, at the most, six 
different eigenvalues. Then we have two cases: 

(i) T has a double and a simple real eigenvalues, and two double and two 
simple complex conjugate eigenvalues iff 2p — 3 3 /3 4 < and a/0, that is 
to say, the Weyl tensor is type IM~ with M 7^ —6. 

(ii) T has three simple and three double real eigenvalue iff 2p — 3 3 /3 4 > 
and /3 7^ 0, that is to say, the Weyl tensor is type IM + . 

Starting from the two cases above with six different eigenvalues we can 
consider three kinds of further degeneration. The first one follows on by 
imposing both conditions, U = tj =/= tk and tj = fj, which is a degeneration 
of both, the IM~ (M ^ -6) and IM+ types. Then, from ([24]) we have 
Tj = f,: = Tj = fj and Tk — U and, consequently, T has four real eigenvalues 
with multiplicities 4, 1, and 4. Moreover, ((23|) implies that the Weyl tensor 
has two equal eigenvalues, p t = pj and, consequently, it is Petrov-Bel type D 
or II. 

Secondly, we can consider a degeneration of the IM~ type by taking the 
three real eigenvalues as equal, t\ = £2 = £3- Then ((24|) implies r\ = t% = 
T3 and T has three different eigenvalues with multiplicity three, a real one 
and a pair of complex conjugates. On the other hand, if we avoid a further 
degeneration, (f2"3")) implies that M = —6. 

Thirdly, we can consider a degeneration of the IM + type by imposing 
Tj = tj. Then ([24|) implies tj — t = Tk = (i,j, k 7^), and T has three real 
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eigenvalues with multiplicities 5, 2 and 2. From (|23|) we have pj = and then 
M = oo. 

Note that, as a consequence of (j!M]) . imposing r, = U to the 7M + case 
leads to vanishing eigenvalues, a more degenerate case that we will consider 
below. All these considerations and proposition |4] lead to the following. 



Proposition 8 (i) A BR tensor T has three real eigenvalues with multiplic- 
ities 4, 1, and 4 iff p? - 4a 6 /? 4 = and 2 2 a 3 = 3 3 /3 2 ^ 0, that is to say, 
the Weyl tensor is type II or D. 

(ii) A BR tensor T has three triple eigenvalues (a real and a pair of 
complex conjugates) iff (3 ^ = a, that is to say, the Weyl tensor is type 
IM [~ 6 1 . 

(Hi) A BR tensor T has three real eigenvalues (with multiplicities 2, 5 
and 2) iff /3 — ^ a, that is to say, the Weyl tensor is type IM°° . 



Finally we can consider a further degeneration from the three cases in 
proposition above. By using (|24|) it is easy to show that in any case we reach 
the highest degeneration: there is a unique vanishing eigenvalue. Then, (|23p 
implies that the Weyl eigenvalues also vanish and, consequently, the Weyl 
tensor is Petrov-Bel type III, N or O. Thus, from proposition 2] we have: 



Proposition 9 A BR tensor T has a sole vanishing eigenvalue with multi- 
plicity 9 iff a = j3 = 0, that is to say, the Weyl tensor is type III, N or 
O. 



It is worth remarking that this proposition implies that the Petrov-Bel 
types ///, N and O can not be distinguished by taking into account the 
eigenvalue multiplicity of the BR tensor T or, equivalently, by analyzing its 
main scalar invariants. We find a similar situation for the Petrov-Bel types 
II and D as a consequence of the first point in proposition [U Nevertheless, 
we will see in the last section that each of these Petrov-Bel types corresponds 
to a different algebraic BR type and that they can be distinguished by using 
the minimal polynomial of T as a TLS-endomorphism. 

Propositions[71[E]and[5]classify and characterize the BR tensor taking into 
account its eigenvalue multiplicity. In order to better visualize the eigenvalue 
degeneration that relates the different classes, we present the following arrow 
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diagram. Note that the four files in the diagram correspond to 9, 6, 3 and 1 
different eigenvalues. 



I, 



{tl, £2, £3, Tl , T2, T3, T\ , Tl , T3} 




IM+ IM 

{tl, tl, t3,Tl,T2,T 3 , Tl, T 2 , T 3 } {t, t, t 3 , T, T, T 3 ,T, T, T3} 





IM°° 

{t,t,0,0,0,T 3 ,0,0,T 3 } 



II, D 

{t,t,t 3 ,T,T,t,T,T,t} {t,t,t, 



IM 1 



T, T, T, 7,7,T 




III, N, O 

{0,0,0,0,0,0,0,0,0} 



6 Obtaining the Weyl tensor from the BR tensor: regular case 



As commented in subsection 12.31 the general problem of obtaining the Weyl 
tensor, up to a duality rotation, in terms of the BR tensor will be solved 
elsewhere [T3] . The specific algorithm strongly depends on the algebraic type 
of the Weyl and BR tensors. Here we present an explicit expression which is 
valid for the Petrov-Bel types D and II and for type / when none of the main 
Weyl scalar invariants vanish. This result is based in a general expression that 
relates the n-powers of the BR and Weyl tensors. On the other hand, this 
expression for n = 2 is used in the last section of this work in showing that 
the minimal polynomial of the BR tensor distinguishes between types II and 
D. 



6.1 Gravitational superenergy tensors of order n 

The tensor T n = W ra o Q for n > 1 is not, generically, a completely sym- 
metric tensor, because W n does not inherit the symmetries of W: it is a 
SD-endomorphism but it is not trace-free (in general) and so the Bianchi 
identity is not satisfied. But its traceless part 



W, 



(n) 



W" 



u> 



(n) 



G- 



oJ 



(n) 



-TrW", 



(25) 
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has the symmetries of the SD Weyl tensor. Thus, it has an associated su- 
perenergy tensor given by: 

T {n) =W(„)oW(„), 

that has all the superenergy properties of the BR tensor T. We will say that 
Ti n \ is the BR superenergy tensor of order n. 

All the expression presented in section [2] involving the Weyl tensor W 
and the BR tensor T will also be valid for W( n ) an d Tt n \ . Specifically, 

T( n ) = f2(n) • i7(„) , 0( n ) = W(„) o G ■ (26) 

From lemma [3] and ([2"5]) the last expression above becomes: 

ft( n ) = fi n - W( n ) r . 

If we use this expression to expand the first expression in (|26[) , we obtain the 
following. 

Proposition 10 The BR superenergy tensor of order n associated with the 
Weyl-like tensor WV n ) defined in (|25p takes the expression: 

T( n ) — T n -W( n ) Q( n )- W( n )J?( n ) - |W(„) | P. (27) 



6.2 A relation between the n-powers of the Weyl and BR tensors 

From the last expression above, we can obtain its self-dual antisymmetric 
part by contracting the indexes (23) with the SD-identity Q. Then, the left 
side vanishes due to the whole symmetry of T/ n ) . On the other hand, the right 
side can be computed taking into account: first, the property of Appendix 
|Pl4, second, the definition (f25|) . and third, that i?( n ) has the same properties 
as Q and, consequently, it satisfies ([12")) and (TIB"]) by replacing W with W(„) • 
Finally, we obtain the following. 

Lemma 7 The power n of a SD Weyl tensor is related to the power n of the 
associated BR-tensor by 



2{T n ) a ^g^ p , 



(TrW")W" --\TrW n \ 2 g 



(28) 

a(3pa 



When Tr W n = 0, then W n = W (n) and T n = T {n) is a completely symmetric 
tensor and, consequently, ([25)1 is an identity. Nevertheless, when Tr W™ does 
not vanish, this proposition provides W™ from T n , up to a duality rotation. 
Indeed, in this case TrT™ = | TrW™| 2 ^ 0. Then, from fl28j we obtain: 

Proposition 11 If T is a BR-tensor and TrT™ ^ 0, the power n of the 
original SD Weyl tensor is given, up to a duality rotation, by 



W n = e ie n 



VTrT' 



^/Tr T n „ 



JSn _ TrW" 

~ |TrW"| ( j 



T(n)appo = {T n )anvp Q^ 'per ( 30 ) 
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6.3 The Weyl tensor from the BR tensor 

Let us suppose now a — TrW 2 ^ and b = TrW 3 ^ 0. Then, for n 
expression (|29| can be written as: 



W 3 



1 



2T ( 



1 



(3) 



\b\ 2 G 



(31) 



Then, if we remove W 3 by using the characteristic equation (|17|) . and we 
take into account that \a\ 2 = TrT 2 , \b\ 2 = TrT 3 , the SD Weyl tensor can be 
recovered from the BR tensor as follows. 



Theorem 1 If T is a BR tensor and Tr T 2 ^ 0, 
tensor can be obtained, up to duality rotation, as 



Tr T 3 ^ 0, then the Weyl 



W 



e W W , 



Wo 



1 



VTr T 2 tr T 3 



4T, 



(3) 



-|TrT 3 |5 



T($) a p P o = (T )a/«/0 £ A 



(32) 



(33) 



This result partially solves the point (III) stated in the introduction. It 
does not apply when af3 = 0, that is, for the classes N, III, IM°° and 
IM'- -6 '. In these degenerate cases the determination of the Weyl tensor in 
terms of the BR tensor requires a different analysis that will be tackled 
elsewhere [13] . 

On the other hand, these results provide an alternative approach to know 
when a traceless completely symmetric rank 4 tensor T is the BR tensor 
associated with a certain Weyl tensor. An answer to this question has been 
given by Bergqvist and Lankien [5] . From theorem Q] an alternative answer 
follows when a and (3 do not vanish. Indeed, if for a given T we calculate the 
SD double 2-form Wo(T) by using expression (f32l) . then T must satisfy the 
equation T = W {T) o W (T). 



7 Complete classification of the Bel-Robinson tensor 

At this point we are ready to complete the classification of the BR tensor by 
considering the minimal polynomial and not only the eigenvalue multiplicity 
(see section [5]). Although a whole analysis of the canonical forms will be 
tackled elsewhere [13] , here we will obtain the necessary results to distinguish 
between types N and ///, and between types II and D, and to intrinsically 
characterize them. 

To accomplish this goal, we start from the conditions on the Weyl tensor 
that distinguish these Petrov-Bel types (see, for example, ref. P35]): 

Lemma 8 Let a = TrW 2 , b = TrW 3 the main Weyl scalar invariants. 
Then, the Weyl tensor is 

i. Petrov-Bel type N iff W 2 = ^ W. 

ii. Petrov-Bel type III iff W 3 = ^ W 2 . 

Hi. Petrov-Bel type D iff a^0,W 2 = A }V- 



^Q- 
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If we apply the property of Appendix [Dl6 to expression (I15p we have 
W™ = if, and only if, T n = 0. Then, considering this property for n = 2,3, 
and as consequence of lemma [8] we can state: 

Proposition 12 Let T be a non null BR tensor. Then, the Weyl tensor is: 
i. Petrov-Bel type N if, and only if, T 2 — 0. 
ii. Petrov-Bel type III if, and only if, T 3 = 7^ T 2 . 

This proposition shows that Petrov-Bel types N and III correspond to 
different algebraic classes of the BR tensor because they have different min- 
imal polynomial. Moreover, it affords a characterization of these Petrov-Bel 
types in terms of the sole BR tensor. The result for type N was given by 
Bergqvist [17] together with the equivalent condition T a p^T" pcrX — 0. This 
requisite states, equivalently, that the superenergy flow vector is a null vector 
for every observer, a condition that was already presented as a characteriza- 
tion of the type N Weyl tensors [16]. Nevertheless, in this work all the condi- 
tions that characterize the different algebraic classes are written as equations 
on the BR tensor as a TLS-endomorphism. 

Now let us deal with a type D Weyl tensor. As the SD Weyl tensor always 
satisfies the characteristic equation (fl7)) . we can compute T 3 = W 3 oW to 
get: 

T 3 = a 2 T + (3 2 T + -{abQ + abQ) . (34) 

On the other hand, from the minimal polynomial of a type D Weyl tensor 

2 

(see point (iii) of lemma [8]) , we can compute T 2 = W 2 o W and we obtain: 

aT 2 = — T+%r+-(aW + abT2), 3 3 /3 2 = 2 2 a 3 . (35) 

3 3^ 6 

Then, we can use ([33]) to remove abQ + abfi in (|3"4"|) and we achieve that 
T satisfies a polynomial of degree 3. Moreover, we know from proposition [8] 
that, for a type D Weyl tensor, the BR tensor has 3 different eigenvalues. 
Thus, this polynomial is the minimal one. The specific calculation leads to 
the following. 

Lemma 9 For a type D Weyl tensor, the associated BR tensor T satisfies 
the minimal polynomial: 

T 3 = aT 2 + ^a 2 T-^a 3 r, a^O. (36) 

We want to prove now that the converse is also true, that is, that condition 
(|3"6"]1 characterize the type D case. Taking the trace of (|3"6"]l we obtain 2 2 a 3 = 
3 3 /3 2 . Moreover, we can use expression (j3"4"]l . valid for an arbitrary BR tensor, 
to remove T 3 in ([36|) . and thus we obtain that T satisfies (f35]) . 

On the other hand, expression (|2"5)) for n — 2 becomes: 

n 2a _ a „ 
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Now, if we compute 7( 2 ) (defined in (j30|) ) putting the expression ((35]) 
for T 2 , and we take into account expressions (JT3J) and (fT5|) and property of 
Appendix |Dj4, we obtain that the SD Weyl tensor satisfies an equation of 
degree 2. Besides, a and, consequently, a does not vanish. Then, the Weyl 
tensor is type D. 

Moreover, if (|55|) does not hold and the main BR scalar invariants satisfy 
the conditions (i) in proposition [HI the Weyl tensor is type 77. These scalar 
conditions can be changed to other equivalent expressions and, finally, we 
arrive to the following. 

Proposition 13 Let T be a BR tensor. Then, the Weyl tensor is: 

i. Petrov-Bel type D if, and only if, T satisfies ((36j) . 

ii. Petrov-Bel type II if, and only if, T satisfies /J, 2 — 4a 6 /5 4 ^ 0. 
2/i = 3 3 /3 4 and it does not satisfy (|36p . 



Again this proposition shows that Petrov-Bel types D and II correspond 
to different algebraic classes of the BR tensor because they have different 
minimal polynomial. 

Propositions [6] [8] [9] [12] and [13] afford a characterization of the differ- 
ent algebraic classes of the BR tensor. These classes correspond to different 
algebraic classes of the Weyl tensor, namely, the Petrov-Bel types and some 
'degenerate' type I cases. All these result can be summarized by the follow- 
ing. 

Theorem 2 For a non vanishing Bel-Robinson tensor T we can distinguish 
nine different algebraic classes. The following arrow diagram offers an al- 
gorithm that distinguishes them. In it, a, (3, /i denote the main BR scalar 
invariants defined in (|20[) . 
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A Products and other formulas involving 2-tensors A and B 

1. Composition as endomorphisms: A- B, 

(A-B)% = A\B% 

2. In general, for arbitrary tensors, • will be used to indicate the contraction of 
adjacent indexes on the tensorial product. 

3. Square and trace as endomorphism 

A 2 = A-A, tr A = A a a 

4. Action on vectors x, y as an endomorphism A(x) and as a bilinear form A(x, y): 

A(x) a — A a gx 13 , A(x,y) = A al 3X a y' 3 

5. Exterior product as double 1-forms: A/\B, 

(A A B) a j3 liv — AcpBpv + Ba^Apv — A au Bp^ — B av App, 



B Products and other formulas involving SD-endomorphisms X 
and y 

1. Every self-dual (SD) symmetric double 2-form X defines a linear map on the 3- 
dimensional SD bivector space. For short, we will say that A" is a SD-endomorphism. 

2. Composition as endomorphisms X o y-. 

{x o y) a i3 pa = — x ^y pa 

3. Square and trace as endomorphism: 

X 2 = XoX, TrX=^X a ^ a/3 

4. Action (on SD bivectors J-, Ti) as an endomorphism X(T\ and as a bilinear 
iormX(T,H): 

X(J-) a f3 = -X a f} J-fiv, X(T,H.) — —Xap^vT Ti 

5. Metric on the space of SD bivectors (SD-identity): 

6. The SD-identity Q satisfies: 

1 — 3 

yafipvy pa = 7Tb/cK/3po-, yafipi^y pa = 7J yaflpa 
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7. The (l,3)-trace of a SD-endomorphism X is proportional to the metric. More 
precisely: 

_ 1 

8. A SD-endomorphism X is traceless iff it satisfies the algebraic Bianchi identity. 
More precisely: 



x»„ = ^Xg aP 



TrX = «• X' , aiip = 4» X aUlltv) =0 

9. The metric volume element r\ is a linear map on the 2-forms space that defines 
the Hodge dual operator. For a real 2-from F and a real symmetric double 
2-form W: 

*F = r)(F), *W = r)oW. 



C Products and other formulas involving TLS-endomorphisms T 
and R. 

1. Every 4-tensor T with the symmetries: 

defines a symmetric linear map on the 9-dimensional space of the traceless 
symmetric tensors. We say that T is a TLS-endomorphism. 

2. Composition as endomorphisms: T • R, 

a/3 rjiap DA^ 



{T.R) al \ a =T a '\ v R»% a 

3. Square and trace as endomorphism: 

T 2 = T»T, Tr T = T ap afl 

4. Action (on trace-less symmetric tensors A, B) as an endomorphism T(A) and 
as a bilinear form T(A,B), 

T{A) a ? = T al T A^, T(A, B) = T aPltv A af3 B^ 

5. Metric on the space of traceless symmetric tensors (TLS-identity): r 

1 1 

r a f^u — - {ga^gfiu + gauges) - jgang^ , r(A) — A 



D Properties involving TLS-endomorphisms defined by 
SD-endomorphisms 

1. For two SD-endomorphisms X and y, we shall define the o— product, X<>y as: 

(x o y) a ij,0v = x a p y f^pvcr 

2. If X and y are SD-endomorphisms, then Q — X O y is a TSL-endomorphism 
such that: 

Q^ fi = \^Qg a ft, TrQ = TrA-TrT 

3. The TLS-identity _T is given in terms of the SD-identity Q by: 

r = g<>g 
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4. The TLS-identity r satisfies: 

3 

-* afivft y pa — "7 y ' a(3pa 

5. If Q = X oQ, X being a traceless SD-endomorphism, then: 

3 

{°Cafi,i/f3 y pa ^ *6 a[pa](3 77 ^\-a0pa 

The last two properties follow from identities of Appendix lBl6. 

6. If A" is a SD-endomorphism, then 

Q = XoX = <S> X = 

7. The SD-endomorphism X is traceless iff Q = X o X is a real fully symmetric 
traceless tensor. 

This property can be easy shown by using the properties 2 above and 8 of 
Appendix [Bl 

8. If V, Z, X and y are SD-endomorphisms, then: 



(VoZ)»(x<>y) = (Vox)o{Zoy) (38) 

This property can be shown by writing in a normalized basis \Ui\ of the SD 
bivectors the four SD-endomorphisms involved, 

and computing the two parts of the expression taking into account that 
g{U u Ui) = Sih V2U t -U 3 =-i e m U k , Q{Ui,Uj) = . 
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